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Weibel instability-mediated collisionless shocks in laser-irradiated dense plasmas: 

Prevailing role of the electrons in the turbulence generation 
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We present a particle-in-cell simulation of the generation of a collisionless turbulent shock in a 
dense plasma driven by an ultra-high-intensity laser pulse. From the linear analysis, we highlight 
the crucial role of the laser-heated and return-current electrons in triggering a strong Weibel-like 
instability, giving rise to a magnetic turbulence able to isotropize the target ions. 


PACS numbers: 

I. INTRODUCTION 

Theoretical and numerical modeling of collisionless 
turbulent shocks is important for understanding various 
high-energy astrophysical environments, where they are 
held responsible for the generation of nonthermal parti¬ 
cles and radiation M. The formation and evolution of 
these structures can now be simulated numerically from 
first principles over significant spatio-temporal scales us¬ 
ing state-of-the-art particle-in-cell (PIG) codes IMil. 
These numerical studies have demonstrated that, for ini¬ 
tially unmagnetized colliding electron-ion flows of high 
enough velocities, the well-known Weibel/filamentation 
instability UMl provides the small-scale magnetic tur¬ 
bulence needed for efficient dissipation of the bulk flow 
energy and Fermi-type acceleration of suprathermal par¬ 
ticles According to simulations, shock forma¬ 

tion proceeds along the following ‘standard’ scenario: 
electron-driven Weibel/filamentation instabilities [illil 
grow and saturate first, leaving the electrons mostly ther- 
malized over the overlap region; for fast enough flows, an 
ion-driven Weibel/filamentation instability subsequently 
develops at larger scales, causing enhanced ion scatter¬ 
ing off amplified magnetic fluctuations; the deflected ions 
then accumulate in the turbulent region, until satisfying 
the shock hydrodynamic jump conditions 

These numerical advances go along with rapid exper¬ 
imental progress towards the generation of collisionless, 
self-magnetized shocks by intense lasers [IB-il- Two 
main configurations are currently investigated to this 
goal. The first one relies upon the interaction of two 
counter-propagating plasma flows generated from the ab¬ 
lation of foil targets by high-energy 0.1 — IMJ), 
nanosecond-duration laser pulses Such flows are of 
relatively low density (n ^ ric, where ric ^ 10^^ cm“^ is 
the critical density of a 1 /im-wavelength laser), temper¬ 
ature (Te ^ IkeV) and drift velocity {vi ^ 1000kms“^). 

An alternative approach, proposed by Fiuza et al [26| 
and further studied in this paper, hinges upon the irradi- 
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FIG. 1: Sketch of the relativistic laser-driven shock formation 
(in the laser piston frame): electromagnetic instabilities are 
triggered in the upstream by counter-streaming electron-ion 
flows. The magnetic fluctuations then isotropize the incoming 
ions, leading to a density jump across the turbulent region 
fulfilling the Rankine-Hugoniot conditions. 


ation of an overdense plasma (n ^ ric) by a relativistic- 
intensity (Jo > lO^^Wcm”^), picosecond-duration laser 
pulse (Fig. [T]). During this interaction, copious amounts 
of electrons are heated to relativistic energies and in¬ 
jected into the target. The ensuing strong charge sepa¬ 
ration and pressure gradients accelerate the surface ions 
to a velocity Vi ^ 2vp, where Vp is the laser-driven pis- 
ton (‘hole boring’) velocity [l^. According to Ref. [26|, 
the various filamentation instabilities triggered in the ion 
beam-plasma region may lead to the formation of a col¬ 
lisionless shock for intense enough lasers and/or dilute 
enough targets. 

In the present paper, we further investigate the phys¬ 
ical scenario of Fiuza et al. [26|. By means of two- 
dimensional PIG simulations, we will mostly focus on 
the laser-specific processes leading to formation of a col¬ 
lisionless turbulent shock, rather than on its subsequent 
propagation. The numerical and physical parameters of 
the simulations will be presented in Sec. [Ill The forma¬ 
tion of a turbulent shock driven by a laser plane wave 
will be evidenced in Sec. imi and shown to fulfill the 
jump conditions of a strong hydrodynamic shock. In 
Sec. la we will demonstrate the dominant role of the 
Weibel/filamentation instability in generating a strong 
electromagnetic turbulence in the upstream region. Us¬ 
ing exact linear theory and approximate non-linear seal- 
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ings, we will assess the respective contributions of the 
plasma species on the instability properties and result¬ 
ing saturated magnetic fluctuations. The prevailing role 
of the laser-accelerated and return current electrons will 
be pointed out, in contrast to the standard astrophysical 
scenario. In Sec. El the ion acceleration and scattering 
around the shock front will be illustrated through a se¬ 
lection of particle trajectories. The late time formation 
of strongly non-linear magnetic vortices will be discussed 
in Sec. ED Shock generation by a focused laser wave will 
be considered in Sec. EID The final section will gather 
our conclusions and perspectives. 


II. NUMERICAL SETUP 

The simulations have been performed using the 
CALDER PIC code, run in 2-D geometry. The laser pulse 
is modeled as an electromagnetic plane wave linearly po¬ 
larized along the ^-axis, with a 1 fim wavelength and an 
intensity Iq = 3.6 x lO^^Wcm”^. This corresponds to 
a normalized field amplitude Aq = eEoIrrieCuoo = 60 
(where ujq is the laser frequency for a Ijnm wavelength). 
The laser intensity is held constant after a linear ramp of 
duration. The wave propagates along the x > 0 
direction and interacts with a fully-ionized, overdense 
plasma slab located at x = dOc/cJo and of maximum elec¬ 
tron density The plasma ions are protons of mass 

mi/rrie = 1836 and charge Zi = 1. The electron density 
is taken to be which corresponds to an elec¬ 

tron plasma frequency of uOpe = \JjmeeQ Iujq. 
These parameters are identical to those considered in 
Ref. [^. The initial electron and ion temperatures 
are = 5keV. A 63c/cc;o (10/im) scale-length den¬ 
sity ramp is added on the front surface to mimic the ef¬ 
fect of the laser pedestal. The simulation grid comprises 
8600 X 1536 cells with mesh sizes Ax = = 0.25c/cc;pe- 

The time step is taken to be At = 0.95Ax/cv/2. Each 
cell contains 50 macro-particles per species, yielding a to¬ 
tal of 1.2 X 10^ macro-particles. To reduce the numerical 
noise, 3rd-order weight factors are used. The bound¬ 
ary conditions are absorbing in x and periodic in y for 
both particles and fields. We have checked that increas¬ 
ing the number of macroparticles to 100 (with a reduced 
box size) or shortening the plasma ramp does not alter 
significantly the results. 


III. GENERAL DESCRIPTION 

Figure shows that, by tujQ = 480, the 

Weib el/filament at ion instability induced by the laser- 
accelerated particles flowing through the bulk plasma has 
given rise to strong magnetic channels parallel to the x- 
axis. From Fig. Hb), we see that these channels are 
associated to pronounced modulations in the ion density 
{dfii/rii ^1). As a result of successive coalescences and 
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FIG. 2: Shock formation in a plasma with Aq = 60 and 
= 50nc. (a) Magnetic field Bz at tuo = 480. (b) Proton 
density nu at tuo = 480. (c) y- averaged proton density vs 
X at successive times. The dashed line shows the initial ion 
profile, (d) y-averaged proton density vs (x,t). 


decreasing anisotropy of the particle momentum distribu¬ 
tions (see below), the field amplitude and filament size 
grow from eBzjrrieUJo 1 and A 3c/uJo (:^ 21c/ujpe) 
to eBzjmeUJo 10 and A 12c/ci;o (— Sbc/oOpe) when 
moving from xujojc 300 to xujojc ^ 200. Closer to 
the target front (150 < xujojc < 200), the filaments ex¬ 
hibit kink-like oscillations, while further growing in am- 
plitnde and size (up to eBzIrrieUJo — 25 and \ujq/c 20) 
. They eventually decay into a compressed turbulent 
layer (further referred to as the downstream) in the inter¬ 
val 140 < xcjo/c < 150. The spatio-temporal evolution 
of this shock-like structure is displayed in Figs. [2jc,d). 
The compression ratio relative to the unperturbed (up¬ 
stream) plasma stabilizes to a value rid/riu — 3 (where u 
and d stand for upstream and downstream, respectively) 
by tujQ 240, and remains constant at later times. The 
right side of the compressed layer propagates at a veloc¬ 
ity '^sh — 0.18c, while its left (irradiated) side is pushed 
by the laser radiation pressure at a velocity Vp 0.12c. 
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(a) tuo = 200 (b) t(jjQ = 240 (c) tuo = 400 

Electrostatic shock phase Electrostatic/magnetic transition Magnetic shock phase 




EIG. 3: Shock formation in a hydrogen plasma with Aq = 60 and = 50nc. Upper panels: y-averaged electromagnetic 
energies normalized to meC^ric. Middle and lower panels: respectively proton and electron x—px phase spaces at (a) tcuo = 200, 
(b) tujQ = 240 and (c) tuo = 400. The dashed contour in the ion phase space at tuo = 400 indicates the phase space region of 
the magnetic vortex studied subsequently (see text). 


These values are consistent with a strong hydrodynamic 
shock induced by the laser-driven motion of the target 
surface at a piston (or ‘hole boring’) velocity 


An'e^rrii 


( 1 ) 


where R denotes the laser reflectivity. This expression 
directly follows from equating the photon and ion mo¬ 
mentum fluxes across the laser-plasma interface (in the 
non-relativistic limit Vp <^1). In the present simulation, 
we have R 0.4, so that Vp 0.12c, which closely agrees 
with the measured value. The theoretical compression 
ratio, rid/riu^ and velocity, Ugh, of a nonrelativistic strong 
hydrodynamic shock are 


'^d Tad T I 

'^u Tad 1 


'^sh 


= Vp 


Tad + 1 
2 


( 2 ) 

( 3 ) 


with Ead the adiabatic index [l9|. In the present 2-D case, 
we have Fad = 2, and hence ndjnu = 3 and Vsh = 0.18, 
which match the simulated values. As first discussed in 
Ref. [26[ , these features bear much resemblance to those 
observed in simulations of Weibel-mediated electron-ion 
shocks in colliding cold flows Q. 


IV. INSTABILITY DEVELOPMENT AND 
SHOCK FORMATION 

The ion and electron x — Px phase spaces around 
the shock front are displayed at successive times in the 
middle and lower panels respectively of Figs. [3fa-c). 


The electrons accelerated in the x > 0 direction have 
a large momentum dispersion {Apx ^ lOOrUec), with 
2cdo-mc)dulations typical of ponderomotive laser acceler¬ 
ation. Later on, the electron distribution broadens with 
time, with maximum px increasing from ^ 400meC at 
tujQ = 200 to ^ SOOrUeC at tooo = 480. 

At tujQ = 200, partial ion reflection occurs off a shock 
front located at xu^jc 105. This gives birth to 
a diluted beam of density ^ O.ln^^, of velocities 
2ush ^ Vx ^ 0.4c, and extending up to xujojc 112. 
In the downstream region, the ions are, on average, ac¬ 
celerated to the piston velocity while exhibiting ve¬ 
locity oscillations. The upper panel of Fig. [3l^a) fur¬ 
ther shows that the ^-averaged electrostatic energy due 
to Ex then slightly dominates the magnetic energy. All 
of these features suggest a transient electrostatic, rather 
than magnetic, collisionless shock driven by the laser ra¬ 
diation pressure [2i-[aii- At tujQ = 240, the reflected and 
upstream ion populations overlap over the space inter¬ 
val 112 < xujqIc < 128. The peak Ex energy, however, 
has decreased by 50% in the vicinity of the shock front. 
Meanwhile, the magnetic energy has increased so that 
it prevails in both the shock foot and downstream re¬ 
gions. Figures SK a,b) plot the 2-D spatial Fourier trans¬ 
forms \Bz{kx^ky)\ and \Ey{kx^ky)\ in the overlap region 
{xujo/c = 120) at tujQ = 240. Both spectra are mostly 
peaked around kxcjuoo ^ 0 and kyc/ujQ ^ 1 — 2, consis¬ 
tently with dominant Weibel/fllamentation modes in the 
shock foot region. The transition into a magnetic shock 
is complete by tujo = 400, at which time the magnetic 
energy exceeds the electric energy by more than one or¬ 
der of magnitude. The shock front has then moved to 
xujqIc 140. Note that the Ex and Ey energies exhibit 
similar profiles in the shock foot region, except at the 
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FIG. 4: Spatial Fourier transforms of Bz (a) and Ey (b) at 
tujQ = 240 and around xujqIc = 120 (in log^^Q scale), (c) 
y-profile of the normalized current density {nVx)lncC of the 
electrons (blue) and protons (red), (d) Normalized filamen- 
tation growth rate F/cco solution of Eq. jnj. The blue line 
corresponds to a multi-waterbag model of the electron and 
proton distributions [Figs. [5][c,d)]; the red line corresponds 
to isotropic Maxwellian protons with Th = 5keV; the black 
line corresponds to immobile protons. 


shock front {xujq/c ^ 140), where the energy becomes 
larger by a factor of ^ 2.5. 

Figure |4l[c) plots transverse lineouts of the electron and 
ion current density, {nvx), at tuo = 240 in the shock foot 
region {xujo/c = 120). The electron current fluctuations 
are seen to largely prevail over the ion current fluctu¬ 
ations. This feature holds at xujojc =115 [Fig. [6l[a)] 
and more generally, on both sides of the shock front un¬ 
til the final simulation time (tcjo = 620). This implies 
that the magnetic filaments, and the resulting shock, are 
mostly driven by the electrons, and not the ions. Simula¬ 
tions performed with /ric) = (60,100), (40,50), 

(20, 50) show essentially similar results. 

This important observation is corroborated by a linear 
stability analysis of the particle momentum distributions 
in the shock foot, displayed in Figs. [5][a,b). As expected 
[ 32 I, the electron distribution is comprised of a diluted, 
high-energy tail (extending up to |p| > lOOrUec) and 
of a denser part carrying moderately-relativistic laser- 
accelerated electrons as well as non-relativistic return 
current electrons. The ion phase space exhibits two 
clearly-separated structures associated to the upstream 
target ions (around 'i; = 0) and the reflected ions (around 
Vx/c 0.35). The latter present a weakly-varying, 
anisotropic (Ap^ {).2mHC ^ SApx) distribution, much 
broader than that of the upstream ions (of temperature 
close to the initial value of 5keV). In order to solve the 
dispersion relation of the Weibel/filamentation instabil¬ 
ity, it is convenient to approximate the measured distri¬ 


FIG. 5: px — Py phase space (in log^o scale) of the electrons 
(a) and protons (b) at xuojc = 240 and tuo = 120. Panels 
(c,d) show the corresponding multi-waterbag fits used for the 
stability analysis (see text). 

butions using a multi-waterbag model [13, [33| 

N 

/(p) = 

i=i 

where aj is the weight of the jth waterbag component. 
The relativistic waterbag distributions are taken to be 
even functions of Py of the form 

Tq(p) = - \p,-P^a\)H{Piy - \py\)5{p,), 

^-Ljx-Ljy 

( 5 ) 

where H is the Heaviside function and Pjd , Pjx , Pjy are 
adjustable parameters. Figures [5][c,d) display the best- 
fitted multi-waterbag approximations (with N = 10^) 
of the electron (c) and ion (d) momentum distributions. 
The multi-waterbag distributions reproduce the momen¬ 
tum fluxes of the original distributions (a,b) to within 
an error of < 10%. The susceptibility tensor associ¬ 
ated to these waterbag distributions is detailed in Ap¬ 
pendix A of Ref. [l3- As shown in Ref. [s^, the 
dispersion relation of the system can be solved for all 
modes co’^(k) G C by generalizing to the relativistic elec¬ 
tromagnetic regime the method introduced by Fried and 
Gould in the non-relativistic electrostatic limit [35|. Us¬ 
ing this procedure, the filamentation growth rate F = "Auj 
is plotted as a function of the transverse wave number 
ky (for kx = 0) in Fig. |4l[d). The maximum value 
Tmax/^o — 0.035 at xujqIc = 120 is obtained at the 
wave number /cmaxc/cJo — 1.3. To assess the contribu¬ 
tion of the ions to the instability, we have replaced their 
multi-waterbag distribution by an isotropic non-drifting 
Maxwell-Jiittner of temperature Th = 5keV. The re¬ 
sulting growth rate curve almost coincides with that ob¬ 
tained with the full ion distribution, which proves that 
the reflected ions are not responsible for the observed fil- 
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FIG. 6: (a) y-profile of the normalized current density 

{nVx)lncC of the electrons (blue) and protons (red) at tcuo = 
240 and xuajc = 115. (b) Normalized filamentation growth 
rate F/oJo solution of Eq. (lAll) . The blue line corresponds 
to a multi-waterbag model of the electron and proton dis¬ 
tributions; the red line corresponds to isotropic Maxwellian 
protons with Th = 5keV; the black line corresponds to im¬ 
mobile protons. 
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(c) xujqIc = 120 





(b) xuojc = 115 



(d) xuqIc = 125 



amentation instability. This feature holds on both sides 
of the shock front and is illustrated in Fig. [6l^b) at 
xuoqIc =115. As will be analyzed later on, because the 
electrons’ mean energy is here comparable to the ions’ (as 
measured in the piston frame), the electrons can induce 
magnetic fluctuations strong enough to scatter the ions 
and entail shock formation. This contrasts with the stan¬ 
dard scenario of Weibel-mediated astrophysical shocks 
j^ . where the unstable two-stream ion distribution is 
considered as the key player in generating magnetic tur¬ 
bulence. Although, in our case, the ion anisotropy is not 
the driving force behind the magnetic buildup, the ther¬ 
mal bulk ions, while stable per se, may significantly en¬ 
hance the electron-driven instability by mitigating space- 
charge effects [Fig. IHd)]. This destabilizing mechanism, 
previously discussed in Refs. [32|, [^, is demonstrated 
by computing the growth rate upon assuming infinite- 
mass ions, and also addressed in [38|. As seen in Fig. 
Ill[d), this yields a maximum growth rate lowered by a 
factor ^ 6, and a dominant wave number down-shifted 
to A^maxC/cJo ^0.7. 

Figures [3[a-d) plot the proton px — Py phase space 
at various locations in the upstream region at time 
tujQ = 240. Both the bulk and reflected proton distri¬ 
butions broaden as one moves closer to the shock front 
as a result of growing magnetic scattering. At the shock 
front {xujo/c = 110), they have merged into a relatively 
isotropized population [Fig. [71[a)]. In Fig. [Sfa) the mag¬ 
netic spectrum \Bz{x,ky)\ in the upstream region is dis¬ 
played and compared to the fastest-growing wave number 
predicted from linear theory using the particle distribu¬ 
tions of Figs. [7l[a-d). Overall, a correct agreement is 
obtained between linear theory and the simulated spec¬ 
trum. 

The spatial profile of the transversely-aver aged mag¬ 
netic field amplitude, (25^)^/^, at tujo = 240 (i.e., the 
approximate formation time of the magnetic shock) is 
plotted in Fig. [Hl[b). The amplitude varies by aboyt two 
orders of magnitude over the upstream region 110 < 
xuoqIc < 140. Note that at this time, the reflected 


FIG. 7: Proton px — Py phase space at tcuo = 240 and various 
locations: xuajc =110 (a), 115 (b), 120 (c) and 125 (d) (in 
log^Q scale). 
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FIG. 8: (a) Magnetic spectrum \Bz{x,ky)\ (in log^Q scale) in 
the upstream region at txo = 240. The triangles plot the 
fastest-growing wave numbers, /cmax, predicted from linear 
theory at the locations of Fig. [7|[a-d). (b) Spatial profile of the 
transversely averaged magnetic field at tuo = 240 

(blue solid line). The black circles plot Eq. (|9]), where feat 
and 7e,eff are measured from the simulation. 


ions extend to xujojc 128 [Fig. [3][b)]. At the shock 
front {xuoq/c 112 ), the magnetic field reaches a value 
2bmeUJole^ comparable to the laser field 

strength. 

In order to analyze this magnetic profile, we have as¬ 
sessed the effectiveness of the various parts of the electron 
distribution in driving the Weibel instability. A similar 
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evaluation was made in Ref. [32|, yet within the simpli¬ 
fying assumption of a purely transverse instability. Here, 
we adopt an alternative approach based on our multi- 
waterbag model. The condition for instability of a multi- 
waterbag system reads A = AB — > 0, where the 

factors given in Appendix B of Ref. 0 , take 

the form of a sum over the waterbag components. For 
instance, we have 


C = J2u^l^P,a/Pfy. (6) 

3 

Let us now define 

Ci = Y,UJ%PjdlP]y, ( 7 ) 

where the index i labels a given waterbag component. 
Likewise, we introduce Bi^ = AiBi — Ci and 
Si = 1 — A^/A. The latter expression then quantifies 
the stabilizing {Si < 0) or destabilizing {Si > 0) in¬ 
fluence of the ith waterbag component. Figures [9l^a,b) 
display SjWj{Y>) using the electron distribution mea¬ 
sured at tujQ = 240, xujqIc = 115 and xu^jc = 120. 
It appears that the electrons mostly responsible for the 
late-stage instability have moderate relativistic energies 
(7 3 at xujqIc = 120 and 7 8 at xbo^jc = 115) and, 

on average, negative x-momenta {px < 0). These parti¬ 
cles are therefore associated with the background return 
current induced by the higher-energy, laser-driven elec¬ 
trons propagating in the x > 0 direction. This feature 
was first pointed out in Ref. [32| under similar laser- 
plasma conditions. The relativistic energies attained by 
the return current electrons stem from the various (elec¬ 
tromagnetic Weibel/filamentation, electrostatic longitu- 
dinal/oblique) beam-plasma instabilities induced in the 
upstream region. 

Let us now confront the simulated magnetic profile 
of Fig. [Hl^b) to simple models of the saturated field ac¬ 
counting for the Weibel-effective part of the electron 
distribution. According to the widely-used transverse- 
trapping model jsol-li^ , saturation occurs when the elec¬ 
tron bounce frequency inside a magnetic filament is equal 
to the maximum growth rate . This yields the magnetic 
field amplitude 

^-^sat / '^ \ f hmax \ ^0 

'^e^O \ (^x / \ ^0 / ^maxC 

where {'y/l3x)efi denotes the average of j/l3x over the 
Weib el-effective electrons. We have typically {'y/Px)ef^ — 
5 — 10. In principle, Fmax is difficult because it is compa¬ 
rable to the laser frequency [H, Q , yielding very short 
exponential growth and saturation time. Using the late¬ 
time Fmax values of Figs. SJ^d) and[6l^b) therefore greatly 
underestimates the saturated field (eHgat/^e^o — 0.01 — 
0.1), except in the far-upstream region, corresponds to 
a more strongly nonlinear regime than is assumed in the 
transverse-trapping model. A second estimate of Hgat 
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FIG. 9: Map of ^SiWi{p) at (a) xujqIc = 115 and (b) 
xujqIc = 120, which measures the local contribution to the 
Weibel instability in the electron px — Py phase space. The 
subpanel zooms in on the most destabilizing electrons. 


may therefore be derived supposing that the Weibel- 
effective electrons are magnetized 1110 . Equating the 
typical filament size, 27r//csat to the Larmor radius of the 
effective electrons, ( 7 )eff/eH, gives the lower limit 


^-^sat / \ ^satO 

- - (7)eff- . 

rrieOJo TTUJo 


(9) 


For a numerical application, the saturated wave num¬ 
ber ksat is extracted from the spectrum of Fig. m (a), 
while the Weib el-effective electrons are defined as that 
part of the electron phase space satisfying SjWj (p) > 
I ma,Xp{^j SjWj). The resulting Hgat values are plotted 
at different locations as black circles in Fig. (H^b), where 
they are found to capture to a good accuracy the simu¬ 
lated magnetic profile. 

In summary, a perturbative analysis of the local plasma 
distribution functions shows that the magnetic turbu¬ 
lence in the upstream region mostly results from the 
interplay between the laser-accelerated and background 
electrons. More precisely, the dominant destabilizing ef¬ 
fect comes from the moderate-energy return current elec¬ 
trons. While the reflected ions have a negligible influ¬ 
ence, the background ions strengthen the instabilit y by 
weakening its inhibiting, electrostatic component 
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(a) 



(b) 



FIG. 10: Typical proton trajectories for various initial x- 
locations: (a) Px{t) and (b) x{t). The color of each curve is 
indexed by Py{t)lmHC. In (b) are also plotted the trajectories 
of the laser-driven piston (black solid line), of the shock front 
(black dashed line) and of the reflected proton front (black 
dotted-dashed line). 


Although the typical wave numbers compare satisfacto¬ 
rily with linear theory, the magnetic field profile in the 
shock-foot region is indicative of the magnetization of the 
Weibel-effective part of the electron distribution. 


V. PROTON TRAJECTORIES 

To gain insight into the dynamics of the upstream pro¬ 
tons, we plot in Figs. flQfa.b) seven typical proton tra¬ 
jectories originating from increasing target 

depths. The color of each trajectory (labeled by the par¬ 
ticle number) is indexed by the instantaneous value of 
the normalized ^-momentum, Py{t)lmHC. In Fig. [TOTb). 
the trajectories of the laser-driven piston, of the shock 
front and of the reflected ion front are also plotted. 

Particles 1 and 2 are rapidly (over the time interval 
190 < tujQ < 220) accelerated forward by the electro¬ 
static field set up by the electron pressure gradient at 
the (then mainly electrostatic) shock front. After reach¬ 
ing the downstream (piston) velocity {vx/c 0.12), they 
remain confined in the downstream region where they 
experience an increasing level of electromagnetic turbu¬ 


lence. The latter makes them oscillate in both Vx and 
Uy, with similar amplitudes Avx ^ Avyjc ^ 0.02. Par¬ 
ticle 3 is more strongly accelerated by the electrostatic 
shock potential: it attains a velocity Vx/c 0.3, which 
corresponds to reflection in the shock frame. During its 
main acceleration phase (220 < tujQ < 300), its ^-velocity 
hardly varies due to a weak magnetic turbulence in the 
shock foot region. Later on, however, the magnetic fluc¬ 
tuations get strong enough to induce velocity variations 
Avyjc ^ 0 . 02 . 

Particles 4-6 exemplify the ion dynamics in the mag¬ 
netic shock regime. Because the upstream electromag¬ 
netic turbulence has grown in amplitude and spatial ex¬ 
tent, they undergo an increasing number of oscillations in 
Vx and Vy while being, on average, accelerated along x by 
the electrostatic field. The effective range of the turbu¬ 
lence can be assessed by noting that particle 7, initially 
located at xbo^jc 165, does not exhibit any significant 
acceleration up to tujo = 500, at which time the shock 
front has moved to xujojc ^ 145 — 150. 


VI. LATE-TIME EVOLUTION: MAGNETIC 
VORTICES 


The late-time evolution of the shock foot region shows 
the formation of partially depleted ion ‘bubbles’. Figures 
IHa- c) zoom in on the bubble located around Xbo^jc = 
180 and yoJo/c = —32, displaying the associated elec¬ 
tron current, transverse electron pressure and magnetic 
field at tuJo = 480. The bubble consists of a mag¬ 
netic dipole sustained by a current of hot electrons flow¬ 
ing in the forward direction, neutralized by thin return- 
current layers at its lower and upper boundaries. The 
mean hot electron current density is je,x ^ —GericC over 
a width ly ^ 7c/uJo , yielding a magneto-static field 
B ^ yolyje,x/‘^ ^ 20meCo’o/e. Thc thickness of the re¬ 
turn current layers is of the order of the relativistic skin 
depth of the upstream plasma ^ ^ O.bc/cJo- 

Similar magnetic structures have been observed in Refs. 
jdol-E^ in the case of lower-density interaction regimes. 
In the present case, the magnetic vortices result from 
the non-linear evolution of the Weibel instability in the 
shock foot region, when the magnetic pressure inside the 
hot electron’s filaments (5^/2/io ^ 200ncmeC^) is strong 
enough to expel the upstream ions. Figure [2fb) shows 
that the magnetic pressure is comparable to, yet lower 
than the hot electron pressure (^ bOOncirieC^). As in 
Refs. the vortex expands at a velocity close to 

the Alfven velocity, Va ^ B/^ 0.1c. The lo¬ 
calized velocity jump indicated by a dashed circle in Fig. 
He) corresponds to a few ions being reflected off a mag¬ 
netic vortex. These magnetic structures develop signifi¬ 
cantly after the ions have become isotropised so that the 
shock formation should not be affected. A few magnetic 
vortices at the early stage of their development are also 
visible in Ref. [26[. Note that these structures, observed 
in our 2D simulation, could be absent in a 3D geometry. 
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(a) Electron current (b) Transverse electron pressure (c) Magnetic field 



FIG. 11: Magnetic vortex: (a) electron current density (nvx) jncC, (b) y-component of electron pressure tensor {nVyPy)IricmeC^ 
and (c) magnetic field eBzIrUeCOQ. 


VII. INFLUENCE OF A FINITE LASER 
SPOT-SIZE 

In order to assess whether the above scenario still holds 
in a realistic setting, we have run an additional simula¬ 
tion with a 16/im (FWHM) Gaussian laser profile. The 
simulation grid has been enlarged up to 16800 x 9216 cells 
with unchanged discretization. The number of macropar¬ 
ticles per mesh has been reduced to 20, so that their total 
number in the simulation is 6.2 x 10^. The other laser 
and plasma parameters are the same as before. 

The simulated ion density and magnetic field at tujo = 
480 are displayed in Figs. fT^a.b) and compared with the 
plane-wave results in Fig. fl^c). Because of the Gaus¬ 
sian shape of the laser beam, the piston velocity peaks 
on the laser axis (^ = 0) and decreases away from it, re¬ 
sulting in a curved laser-plasma interface. However, the 
proton density of Fig. [T^ a) shows a relatively flat trans¬ 
verse profile (independent of ^), close to the laser axis 
{\yuJo/c\ < 20). As further shown by the density line- 
outs of Fig. [T^c). the plasma is compressed by a factor 
^ 3, consistently with the Rankine-Hugoniot conditions. 
This compressed region exhibits about five magnetic fila¬ 
ments, which seems to be sufficient to form the shock-like 
structure. Lineouts of the proton density are plotted at 
different times in Fig. [T^ ch and compared to those from 
the plane wave simulation. At tuJo = 240, the two con¬ 
figurations yield similar profiles. At tuJo = 480, the lo¬ 
cation of the laser-piston and the compression factor are 
very close in the two cases, yet the downstream region is 
thinner (by ^ 50%) in the focused configuration. This 
can be explained by the curvature of the laser-plasma in¬ 
terface [Fig. [2l[a)], which allows the ions and electrons 
of the downstream to leak off the y = 0 axis. To en¬ 
sure a stable front shock over a duration At, the laser 
spot size, should be larger than VpAt, otherwise, the 
downstream will be significantly distorted by the Gaus¬ 
sian piston, which may hamper the shock propagation. 
In the case of a D = 16/im-focal spot, the shock should 
be stable up to At < D/vp 800co’^^. The drop in 
the downstream density observed at tujQ = 620 (corre¬ 
sponding to effective interaction time of 520 cl;^^) sup- 


(a) uhItIc (b) eBzIrUeOJo 




FIG. 12: (a) Proton density, (b) magnetic field at tuo = 480 
and (c) ion density profiles averaged over the transverse direc¬ 
tion in the region \yuJo/c\ < 10 (red solid line) for a Weibel- 
mediated collisionless shock induced by a laser of intensity 
Iq ~ 1.8 10^^W.cm“^ and 16/im spot size. The dashed black 
lines plot the //-averaged density from the plane wave simula¬ 
tion. 


ports qualitatively this estimate. Note that making use 
of a super-Gaussian laser profile may yield a more stable 
shock front. 
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VIII. CONCLUSIONS 


tions yields the well-known dispersion relation 


In contrast to the currently explored experimental 
setup [ 2 ^, [ 23 , Eil, [ 5 ^ which requires 10 —lOOkJ-class laser 
facilities, we have found that only 1 — 2kJ of intense- 
enough laser pulses should be able to drive Weibel- 
mediated shocks, which confirms the results of Ref. [26|. 
Their formation proceeds as follows. Rapidly after the 
start of the irradiation, an electrostatic shock forms and 
propagates. At later times, because the hot electrons 
have typical energy comparable to that of the ions (in 
the piston frame), they are able to trigger a strong mag¬ 
netic turbulence in the shock front region. The crit¬ 
ical role of the return-current electrons in driving the 
Weib el-filament at ion instability has been shown, making 
use of a multi-waterbag decornposition scheme. Unlike 
the usual framework [Sbl, EH, [5l|, the ion heating and 
isotropization result from an electron-driven instability 
(rather than ion-driven), at least over the relatively short 
time-scales considered in our simulation. During shock 
propagation on longer time-scales, the ions could play a 
significant role in triggering the upstream instability. We 
have also shown that shock formation is accessible to a 
broad enough focused laser wave. However, the shock is 
found to decay away over time scales > Djvp 

Several points have yet to be clarified to ensure the ex¬ 
perimental feasibility of Weibel-mediated laser-induced 
collisionless shocks. The influence of collisions, ioniza¬ 
tion, radiative losses as well as 3D effects could be im¬ 
portant. One of the most critical issue is linked to the 
diagnostics of such experiments. A major drawback of 
such shocks is that they develop in dense plasmas over 
very short space (^ /im) and time scales (lOOfs), which 
greatly complicates their experimental characterization 
unless one disposes of an intense, short-duration x-ray 
probe as provided by a free-electron-laser [52|. 


{ujSxx — sin^ 0){ujSyy — cos^ 0) 

— {uJSyz — k‘^ cos OsiiiOY = 0 . (Al) 

The dielectric tensor reads: 


s 

where Y is the relativistic susceptibility tensor given in 
Refs. [l3,0lll| for waterbag and Maxwell-Jiittner dis¬ 
tribution functions. All the tensor elements can easily be 
recast as a function of the phase speed and the angle 0. 
Equation dSB) gives 

ak^ + bk^ -1- c = 0 , (A3) 


with 

a = (/^^ — sin^ ^)(/50 — cos^ 0) — cos^ 0 sin^ 0 , (A4) 

b = {sin^e- 0l)'^u)l^xyy{cos^ 0 - Pl)^u)l^Xzz 

s s 

+ 2 cos 0 sin 0 E ^IsXyz , (AS) 

C= C^<^lsXvv)C^<^lsXzz) - C^^lsXyz? , (A6) 

s s s 


In Eqs. (jA4p . (jA5p and (IA6p . the subscripts have been 
omitted on y for the sake of clarity. The wavevector then 
verifies 


7 2 ^ -6(/30) -h 

MM 

2 _ - 

2a{M 


(A7) 

(A8) 
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Appendix A: Dielectric tensors 

Introducing 6>, the angle between the wavevector and 
the x-axis, the linearization of the Vlasov-Maxwell equa- 


with A = \/h‘^ — 4ac. This formulation, in which /c^(> 0) 
is a function of only, lends itself to the efficient numer¬ 
ical scheme introduced by Pried and Gould [35| in a non- 
relativistic framework and generalized recently [sj • This 
scheme consists, first, in determining the locus of the ze¬ 
roes of This can be readily performed by means 

of a contour plot in a finely discretized portion of the 
complex 13(1) plane. Then, we retain those zeroes fulfill¬ 
ing ^G{l3(p) > 0 and identify k = -s/WgIM)' Depending 
on the Pcj) domain considered, this method allows us to 
simultaneously solve for a set of discrete electromagnetic 
solutions uj{k^ 0). 
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